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AUSTRACT 


A computer  program  has  been  written  in  Fortran  language  tliat  plots 
fictitious  latitudes  and  longitudes  about  an  arbitrary  polar  system  onto 
standard  Mercator  projections.  The  matitematical  development  of  tliis 
program  Is  presented  along  with  a complete  program  description  and 
program  listing.  Several  example  plots  generated  by  the  program  are 
included  to  demonstrate  its  option  cliaracterlstics . 
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INTRODUCTION 

This  report  describes  a computer  program  that  plots 
latitudes  and  longitudes  about  any  fictitious  polar  system 
onto  standard  Mercator  grids  of  geographic  latitudes  and  lon- 
gitudes. The  principle  of  this  program  is  shown  schematically 
In  Figures  1 and  2.  Figure  lA  shows  the  geographic  latitude 
and  longitude  system  for  one- ha  If  of  the  earth,  0*  to  180° 
longitude.  Figure  IB  1s  the  standard  Mercator  projection  at 
10°  Intervals  of  latitudes  and  longitudes  for  70°N  to  70°S 
latitude,  and  0*  to  180*  longitude.  In  Figure  2A,  one-half  of 
the  earth  is  shown  schematically  again  with  geographic  latitudes 
and  longitudes,  only  this  time  with  tlie  addition  of  a second 
latitude-longitude  system  for  a second  polar  pair,  Pj  (60°N, 

180°)  and  (60°S,  0°).  In  Figure  2B,  the  fictitious  polar 
system  of  latitudes  and  longitudes  Is  shown  projected  onto  a 
standard  Mercator  grid. 

The  mathematical  development  of  the  program,  presented  In 
the  following  section,  should  allow  the  reader  to  adapt  the 
program  to  other  computer  languages,  such  as  PL/1,  Assembly, 
or  Basic.  The  program  description,  which  includes  detailed 
Instructions  for  the  use  of  job  control  cards,  allows  the  program 
to  be  utilized  by  workers  with  limited  experience  on  an  IBM/360 
computer.  A complete  program  listing  Is  Included  In  the  final 
section  of  the  report,  and  examples  of  several  plots  generated 
with  the  program  are  presented  in  the  Appendix. 


Fig.  1.  A.  Global  scheaatic  of  standard  longitudes  for  one-half  of  the  earth  (0»-180«  longitude) 
B.  Mercator  projection  for  70'N-70“S  latitude  and  0“-180*  longitude. 
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A. 


Although  this  program  was  developed  for  utilization  In 
geotectonlc  studies  seeking  to  describe  the  surface  kinematics 
of  the  earth's  crust  (for  an  example  see  Figure  15  of 
Handschumacher , 1976),  its  applicability  to  numerous  other 
scientific  investigations  is  apparent.  The  authors  are  inter- 
ested in  learning  of  such  application  anil  would  appreciate 
receiving  Information  on  how  this  program  was  used  in  other 
a reas  of  study  . 


c>u 


1 


I 


I 


MATHEMATICAL  DEVELOPMENT 

Consider  tliat  the  axis  of  the  cylinder  tangent  to  the 
sphere  of  radius  R points  in  the  direction  of  the  unit  vector 

K = i SIN  COS  (|)  + j SIN  e SIN  (|)  + k COS  0 , (1) 

'v'V.  p P % P 0"“  P 

where  0 Is  the  polar  angle  of  the  axis  and  (J)  is  the  azimuth 
P p 

angle.  Here  the  z-axis  is  the  earth's  axis  passing  through 
the  poles. 


Fig.  3. 

We  may  find  two  other  vectors  perpendicular  to  and  to  one 

another  by  using  an  analogy  with  the  spherical  coordinate 

mutually  perpendicular  vectors.  If  K corresponds  to  e^  and 

% 

to  e and  I to  e , we  have 
a.cj)  a/  a,0 


6. 


1 = 1 COS  0 COS  4)  + j cos  0 SIN  * - k SIN  0 

^ P p p P 'X.  p 


.I  = -iSINd)  + i COS  (b 
% O,  P p 


(2) 


K = 1 SIN  0 COS  4i  + j SIN  0 SIN  d>  + k COS  0 
P pK,  P P%  P 


I,  J,  and  K establish  another  rectancular  coordinate  system 

'h  % 'X,  ^ 

X,  Y,  and  Z as  shown  In  Figure  3. 

Every  point  on  the  cylinder  of  radius  R with  axis  Z has 
position  vector 


r = R COS  t I + R SIN  <I>  J + Z K , (3) 

<V,C  "X/  'V  % 

where  $ ^Is  the  azimuth  angle  in  the  X,  Y,  Z system.  The  radial 
line  with  parameter  t with  azimuth  angle  4>  the  same  as  equation 
(3)  is 


r = (SIN  0 COS  4)  I + SIN  0 SIN  4>  J + COS  0 K)  t . (4) 

% 

Equations  (3)  and  (4)  Intersect  when 
t SIN  0 COS  4>  - R COS  4> 

t SIN  0 SIN  4>  - R SIN  4>  (4a) 

t COS  0 - Z . 


i'q  nation  (4a)  gives  us 


t SIN  G = R 

t COS  G ' Z 

or  COT  G = Z/R 

t = R/STN  0 

Thus  from  equation  (3)  and  (4)  equivalently,  we 
equation  (5) 

r„  = R/SIN  0 (SIN  0 COS  4>  I + SIN  G SIN 

'C"  % 

•=  R COS  <tl+RSIN<I>J  + KZ 

'll  % 

COT  0 » Z/R  . 

Now  if  we  consider  the  cylinder  of  radius  R 
is  the  z-axls  and  the  radial  line  with  parameter 
similar  to  the  above  in  the  x,  y,  z system 

r * (SIN  e COS  * i + SIN  6 SIN  ()>  1 + COS 
<v.m  'x,  'll 

r = R COS  (J)  1 + R SIN  <|)  J + z k 

<\,C  'll  h 

SOj  similar  to  equation  (6),  the  position  vector  i 


obta in  us Ing 
J + COS  0 k) 

'V  'll 


whose  axis 
S , ..we  have 

0 k)S 

'll 

f inter- 


section is: 


8. 


r * R/SIN  0 (SIN  9 COS  <)>  1 + SIN  0 


SIN  (j)  J + COS  9 


k) 


- R COS  (J»  1 + R SIN  (|»  J + k z 

'Vj  'll  % 


(8) 


COT  e - z/R 


Now  la  parallel  to  r since  we  are  trying  to  find  a 
correspondence  with  points  on  the  sphere  to  its  prolectlon  on 
the  cylinders  by  radial  lines.  Since  r„  Is  parallel  to  r , 

'V,n  fY,m 

we  must  have 


In 


C r 


(where  C Is  a constant). 


(9) 


Thus, 


Cr  . So  from  equation  (6)  and  equation  (8) 

(Vira  * 


R . P R p . SIN  9 

SIN  0 ^ SIN  0 or  L siN  0 * 

Therefore,  equation  (9)  may  be  written 

SIN  0 r„  - SIN  9 r . (10) 

'\,n  <\/m 

On  Introducing  the  second  equations  of  (6)  and  (8)  In  equation 


( 10) , we  obtain 


I 


9. 


SIN  0 (R  COS  4>  I + R SIN  4>  .1  + K Z) 

% '\y  a, 


SIN  ft  (R  COS  (j)I  + RSINit>i+k7. ) 


If  equation  (11)  Is  dotted  In  turn  by  1,  j,  k,  with  the 

% <\j  % 

use  of  equation  (2),  we  obtain  the  following: 


(11) 


R SIN  9 COS  (t,  » SIN  0 (R  COS  <I>  COS  6 COS  * 

P P 


- R SIN  4)  SIN  d)  + Z SIN  0 COS  d>  ) 
P P P 


R SIN  ft  SIN  A * SIN  0 (R  COS  4'  COS  0 SIN  A 

P P 


(12) 


+ R SIN  4>  COS  (j»  + Z SIN  0 SIN  (p  ) 

P P P 


z SIN  0 - SIN  0 (-R  COS  4>  SIN  0 + 0 + Z COS  0 ) 

P P 


If  the  first  two  equations  above  are  squared  and  added 
together,  we  obtain 


SIN  0 - SIN  0 (COS^  4>  COS^  0p  + SIN^  ♦ + (Z/R)^  SIN“'  9^ 


+ 2(Z/R)  COS  A SIN  0p  COS  0^) 


1/2 


(13) 


Thus  from  equation  (12)  we  obtain 


10. 


cos  cos  0p  cos  4>p  - SIN  4>  SIN  (|)p  + (Z/R)  SIN  0p  SIN 

COS  *)*  “ / 2 2 2 2 

/ COS  'I' COS  0 + SIN  4>+  (Z/R)  SIN  0 + 2(Z/R)  COS 'I' SIN  0 COS  0 

/ P P p p 


SIN<t> 


COS  <1*  COS  0 SIN  <l>  + SIN  COS  + (Z/R)  SIN  0 SIN  ^ 

P P P P P 

cost'll  COS^  0 + SIN^4>+  (Z/R)^  SIN^  0 + 2(Z/R)  COS ♦ SIN  0 COS  9 

P P P P 


(14) 


z/R 


- COS  $ SIN  0 + (Z/R)  COS  e 

P P 


COS^  <I>COS^  0 + SIN^  4>  + (Z/R)^  SIN^  0 + 2(Z/R)  COStSIN  0 COS  0 

P P P P 


Equarlon  (14)  enables  the  (4),  z)  to  be  expressed  In  terms  of  (♦,  Z) , 
Let  us  make  the  following  conventions: 


NP 


where  0 - latitude  of  fictitious  pole  referring  to  True  North 

P 

Pole 

- longitude  of  the  fictitious  pole  referring  to  True 
North  Pole 

4>  - 0 passes  through  Greenwich 


11. 


9 =0  North  Pole 

P 

9 = TT  South  Pole 

P 

is  measured  from  west  to  east. 

Thus  53®W  has  for  <J)  , <j>  = 360"  - 53*  . 

P P 

Also  53"N  corresponds  to  6^  ■ 90*  - 53*  . 

Notice  that  if  the  fictitious  pole  is  south  of  the  equator 
then  53*S  would  correspond  to  9^  ■ 90*  + 53*. 

The  simplest  way  to  generate  the  map  of  fictitious  longi- 
tude and  latitude  curves  on  the  Mercator  map  is  to  use  equation 
(14)  . 

Define 


A = COS  4)  COS  0 COS  d)  - SIN  4>  SIN  d»  + COT  0 SIN  0 COS  4 

p '♦'p  '^p  p '<'p 


B S COS -t  COS  0 SIN  A + SIN4>COS  d>  +COT0SIN0  SIN  A 
P ^P  P ^P 


(15) 


R = / COS^  4i  COS^  0 + SIN^  4»  + COT^  0 SIN^  B + 2 COT  0 COS  « SIN  0 COS  0 

/ P P P P 


2 2 


where  we  have  put  Z/R  ■ COT  0. 

To  generate  the  maps  of  the  fictitious  longitudes  put 

<I>  “ <ti  into  (15)  where  4>  is  varied  over  the  fictitious  merl- 
P P 

d 1 a n s . 

Then  from  equation  (14)  we  have 


12. 


(J)p  = TAN“^(B/A)  or  SIN~^(B/(A^  + B^) 

- COS  $ SIN  0 f COT  0 COS  0 
COT  0 - z/R  P ■ P P (16) 

/ 

and  vary  © through  0 < 0 < tr  to  generate  the  map.  Equation 
(16)  is  the  parametric  equation  for  the  map  with  parameter 
being  0 . 

To  generate  the  maps  of  the  fictitious  latitudes  we  put 
0 ■ 0p  into  equation  (15)  above  and  let  range  through 
0 < ♦ < 2 Tf  to  generate  the  maps  corresponding  to  each  0^  . 

_ 1 g 

Note  that  if  Tan  -r  is  used  we  must  put 

A 

♦ - TAN"^  (B/A)  + (1  - H(A)]  * TT  + h(A)  * [1  - H(B)]  * 2 


whe  re 


H(X)  - 0 If  X < 0 
- 1 If  X > 0 

However,  If  the  above  convention  is  not  used,  the  follow 
Ing  can  be  considered: 

(ji  - TAN"^(B/A)  < 0 WEST 

> 0 


<P  - TAN“^(B/A) 


EAST 


13. 


PROGRAM  DESCRIPTION 

The  program  listed  in  this  report  is  written  in  Fortran  IV 
at  a sufficiently  generalized  level  to  be  compatible  with  or 
easily  adapted  to  most  computers  and  plotters.  It  has  been 
compiled  and  executed  successfully  on  IBM  370  and  IBM  360  com- 
puter systems,  and  the  plotting  has  been  tested  on  Calcomp  and 
XYNETICS  plotters  as  well.  The  program  is  subdivided  into  a 
Main  program,  and  three  subroutines,  FICT,  EDGEPT,  DEGREE,  and 
also  a function  program  YMER. 

The  function  of  the  Main  program  is  to  enter  the  input  data. 
It  computes  the  required  constants,  converts  the  coordinates, 
and  checks  the  map  for  the  size  in  Y-axis.  It  calls  in  sub- 
routine DEGREE  to  develop  lat 1 1 ude- long  1 1 ud e annotation  at  the 
left  and  the  bottom  for  the  Mercator  map.  The  arrangement  of 
data  cards  and  various  options  are  discussed  in  the  section  on 
data  input  preparation.  Once  the  Mercator  map  has  been  com- 
pleted, it  calls  the  subroutine  FICT  to  generate  all  the  ficti- 
tious curves  on  that  particular  Mercator  map,  and  last,  it 
either  looks  for  another  sot  of  data  for  Pole  projection  or 
terminates  its  operation. 

Subroutine  EDGEPT  is  called  by  the  subroutine  FICT  to 
extrapolate  the  fictitious  curve  to  the  edge  of  the  Mercator 
map  if  a data  point  is  being  found  outside  the  map. 

Function  program  YMER  is  called  by  the  Main  program  and 
subroutine  FICT  to  calculate  the  displacement  from  the  equator 
to  transform  the  point  found  on  a sphere  Into  a point  as  if  on 


an  ellipsoid. 


lA. 


Dimensions 

The  program  Is  dimensioned  so  that  fictitious  longitude 
and  fictitious  latitude  curves  are  obtainable  for  up  to  160 
parameters.  This  number  Is  arbitrary;  however,  it  may  be  made 
larger  by  increasing  the  dime'  :<lons  for  tlie  appropriate  program 
variable,  FPL ON.  No  other  changes  in  the  program  are  necessary. 

If  one  wants  to  process  additional  plots  In  one  execution 
time,  another  set  of  Input  cards  is  necessary  and  should  be  In- 
cluded, as  shown  below,  after  the  comment  card  of  the  previous 
set  of  data . 


//MERPLOT  JOB  (9999,  IM,  3KI,  3KL) , PETER 
//  EXEC  FORTCLG 
//SYSIN  DD  * 

FORTRAN  SOURCE  PROGRAM 

//GO. PLOTTAPE  DD  UNIT- ( 7TRK ,, DEFER)  , LABEL- (, NL)  , 
//  VOL-SER-PLOT,DISP-(NEW,KEEP) 

//GO. SYSIN  DD  * 

AO.OON  110. OOW  10  10 


50N 

SOS  70W 

130E  .063 

10 

1 

1 

POLE 

PROJECTION 

AON  HOW 

/* 


15. 


Minor  modtflc.Tflons  to  the  prof»rani  are  necessary  for  a 
successful  compilation  and  e x e e ii  t i j) n on  the  o t li e r computers. 
Statements  like  CAI.I.  IM.OTS  , CAM.  PLOT,  CALL  SYMBOL,  have  to  be 
changed  If  the  plotter  used  (s  not  the  Calcomp  or  XYNF.TICS 
plotter.  The  size  of  Y-maximum  can  be  Increased  more  than  the 
lO-inch  limit  If  the  allowal)le  size  on  Y-axls  Is  bigger  than 
that  of  the  11-lnch  C.alcomp  plotter. 


Input  Prep a r a 1 1 on  and  Option s 


CARD  1 

READS  IN  LOCATION  OF  THE  POLE  AND  THE 

FICTITIOUS  GRID  SPACINC  ( DECREE/ INTERVAL) 

cOL. 

1-5 

LATITUDE  OF  POLE 

F5.  2 

COL. 

6 

NORTH  OR  SOUTH  (N  OR  S) 

A1 

COL. 

8-13 

LONGITUDE  OF  POLE 

F6 . 2 

COL. 

lA 

EAST  OR  WEST  (E  OR  W) 

A1 

COL. 

15-17 

GRID  SPACING  IN  FICTITIOUS 

LONGITUDE  CURVES 

F3.0 

COL. 

18-20 

GRID  SPACING  IN  FICTITIOUS  LATITUDES 

F3.0 

CARD  2 

READS  IN  THE  SIZE  OF  THE  MAP 

COL. 

1-5 

TOP  OF  THE  MERCATOR  MAP 

F5 . 0 

COL. 

6 

NORTH  OR  SOUTH  (N  OR  S) 

A1 

COL. 

8-12 

BOTTOM  OF  THE  MERCATOR  MAP 

F5.0 

COL. 

13 

NORTH  OR  SOUTH  (N  OR  S) 

A1 

COL. 

15-19 

RIGHT  LONGITUDE  OF  THE  MERCATOR  MAP 

F5. 0 

16 


COL.  21  EAST  OR  WEST  (E  OR  W)  A1 

COL.  22-26  LEFT  LONGITUDE  OF  THE  MERCATOR  MAP  F5.0 

COL.  27  EAST  OR  WEST  INFORMATION  (E  OR  W)  A1 

COL.  29-34  MERCATOR  SCALE  (INCH/DEGREE)  F6.3 

COL.  36-41  MERCATOR  GRID  SPACING  (DEGREE/ INTERVAL)  F6 . 3 

COL.  55  IF-1  FICTITIOUS  LONGITUDES  ARE 

DRAWN  ON  THE  MERCATOR  MAP  II 

COL.  59  IF-1  FICTITIOUS  LATITUDES  ARE 

DRAWN  ON  THE  MAP  II 

CARD  3 

COL.  1-64  TITLE  OF  THE  PLOT  16A4 


n n r. 


17. 


F’ROCRAM  LISTING 


c 

c 

c 

c 


TMi*-.  (■i  i'(,t  *M  in  Tc  i‘i  I T ricTiTimin  mh'Ioians  anh  paraiitl’ 

RM.|i|*r,  MbPl-ATl’U.  »»»  NllV.  II  |Q74,  ••• 

HAWAII  insitrotr  (If  r.f^ni’HYSi cs 


ON 


r dmkiin  amt  . kmm  . vmt  , y Mn  . V rm  , vmiT  .ssr  ,A  P I r.f. 
INK  (.n  A,it,(.,(i,N/»N*/,r/*F*/ 

ij  I Ml  nm  ON  III  If  I ( I oni>  I , T I T I r ( I M 

c 

c • * • • • I N I T I M I / r T HI  0|  (IT 

V 

(_»»»»•(_  A|  I ' ■(  I I ( 1 , I . 

C»»*«*(All  I'l  II  T ',  ( HI  II  I 

c 


11  I . I M T 1 AL  <•  A| 

. A OOP  I Is  I M T I A1 


nip 

F Al  L 


■ T I T I 


xvNi  T irs  fM  ri T Fo 

row  fALtOMP  (»|  rjTTf  P 


I 4 
M 


10 


t Ai  L I’l  m s < mil  r . aoou  i 

Ai'ni  - n . 

V A » ip  . 

(•fAH  (n,A.iNi)=7i  ri  A T , NS  ,r  lon.ew  .nru  CN.  r.CL  AT 

fill- MAI  (I  Ai  , 1 *.rf>,.i,Ai  ,r  1.0, r 5.p> 

Pf  AH  ( •,  . 1 (1  ,1  NO  1 7 I V T . A , VH  . M . *R  .C  . XL  . n .SCAl.f  ,I.P  ID  , 1ST 
Flll,MAT(4(rS.().Al,|«|,?<rF,  .I.lXI.IOX.TlII 
p f AIK  S , T ,1  NO  -T  1 T I Tl  r 
r UP  M A T 1 I 7.A  A I 

PRINT  II  ,y  T . A . vM ,n , xP ,r , XL .n. seal  r .SRI  0 

f 111  MAT  ( • I • . • I'l  OT  limits:  • .4 1 rr . j.a  1 , 2x  )// ,ax  , • scaie-'.ft. 

I.  • II  or  S/l'Ti.l  I r • // , TX  . •i.Rlo  SPACING  **.rT.T,»  or  f.REr  S/ I NTFRVAl 
I’l-  I II  I 1 , I ST 

r u ( • M A 1 ( • o • . ■ « , t I 2 X A s f<  ''  •/ 
f.  Tx,M||,?«|/ 

. pd.TST  sA  •)',  so  ‘,r  SH  SQ*i 

PRINT  I 2 , II  Tl  f 
roPMAIl///SX,lf.AA//l 


■ fiPINT  MAP  Tiur  AND  SFT  THf  PIN  TO  ORIGIN. 


C A I L syMIlut  ( «I>(IT  *2  ,S  , V Ax  . . I A ,T  I Tl.r  ,0  .,f.A  I 

call  I'l  dm  xhhi  ♦ ? , 0 . V ax  AO  . rs.  1 

SSC  =SCAl ( /hO, 

I r ( A.  ni. N.  and,  n.  Fo.N I go  TO  t 
IF  I A.NE.N.AMn.M.NF ,N|  GO  TO  2 
YM  T S y T»  fiO 

YMfi  =-  yn»(jO 

GO  Tl)  1 
YMT  =Y  T*G0 
YMfIsYHAf.O 
GO  TO  3 
YMT  =-  YT4A.0 

YMMi- yn«  60 

I f I C.  I O.  t . AND,  O,  CO.e  I GO  TO  A 
It  I C . NF  , F . ANO.  O.  NF  ,r  I GO  TO  5 
XMT  ■=!  0000.  -XI  •GO 
XMM=XL«ftC- 10000 . 

GO  TU  t 

XMT -XP4G0- 10000 . 

xMii  :xL  AGO-  inooo . 

GO  T O 6 

XMT il onoo. -XPA60 
XMM  = 1 onoo. -XL*  GO 
YM0T= YMFRI YMHI 
YTrP3<YMrR|VMT)-YHOT)*  SSC 

‘CHECK  Tlir  OVFRALL  SI/E  OF  THF  MAP.Y-MAX,  IS  10  INCH, 

> THf  RF  IS  NO  X-MAX.  ON  THF  CALCOMP  PLOTTTR  T E C HN  I C Al  L Y , 


IF  lYTOP.Gf  .10.  I GO  TO  33 
XR  If,M  = ( XMT-XMh)ASSC 
PRINT  17,  XRIGH.YTOP 
17  FORMATI/SX,'  MAP  51/E  IN  INCHES* 
t*  X : • ,r  7. 1 ,sx  ,*  Y : • ,r  7.3////I 
NX  Ax  I XMT-XMM I /(  GO. *GP I Ol *1  «S 
XMx(  XMT-XMITI/INXA-I  I 

no  21  iai.nxa 

XPOTXXMA  n -1  I ASSC 

Xf  *|7XH()»  ( I -I  I AGH  I DA60, 

ir I Mool I .21 .ro.o I go  to  77 


C AAAA  AORA  tt 

c 


VFRTICAL  GPIOS  ANO  ANNOTATE  IN  DEC.REES. 


77 


CAIL  DEGM' r ( XPOT  ,XF«I.  I 
CAl  L PI  (IT  I XPilT  .0.  . Tl 
CALL  PI  t;T  I XPOT  , Y TOI’.  2 » 
GO  Til  21 

CAIL  PI  OT  I XPOT  ,Y TOP. 3 I 
CALL  *’l  OT  I XPOT  ,0  . » 


**  * »OOOiK  n 1 r, 

• oomr  020 

•nnrnfn3c 

•noOOPPAO 

« no  OOP  OGO 

Aooonpnr  o 

aaw.00000070 

oonoccpo 
ooooooso 
000001 00 
ooror 1 1 0 
noooo I 20 
oonno I 30 
prooo 1 AO 

00000 1 FO 

ocnoo 1 FO 

000001 70 
000001 »0 
000001 00 
0000O200 
0000021 0 
00000220 
00COO2  30 
000002AO 
OOOOC2FO 
000002FO 
POOCOi'  70 
3.  ooooorno 
•221000002<5C 
00000300 
O0COO.3IO 
00000320 
POOOO  3 30 
0O0O03A0 
00000X50 
00000360 
00000370 
OOOOC  3E0 
00000  xi;o 
OOOOOAOO 
oooor Al 0 
OOOOOA20 
OOOOPA  30 
OOOOOA40 
00000450 
00000460 
000004  70 
O0C004PO 
00000490 
OP0OO50P 
COOCOSIO 
00000520 
00000*  *0 
00000540 
00000550 
00000560 
00000570 
00000560 
00000590 
00000600 
0000061 0 
00000620 
OC000630 
00000640 
00000650 
00000660 
00000670 
00000660 
OOOOC690 
00000700 
0000071C 
00000720 
00000730 
00000740 
00000750 
00000760 
00000770 
OOC00760 
00000790 

ooooonoo 

ooooonto 

00000820 

00000030 

OOOOOPAO 

00000850 

00000860 

00000870 

00000880 


nnon  nrinnrionn  onn  oor.n  -rir. 


IH 


c»«»»»ir  ■'.>■  Ai  r T np  '.ma(i  4^.>  (.f  in  ; Tim  n o^r  fvrr  v nTiin: 

t . • • . I 1 n (|,  I ,,|  I p j ANMin  A Tf  II  Til  AVI  I I II  | II V I I II  | | | T Y 

I ' * •'  **  • • I ’ . o . oi  . and . (,i  1 1) ,1  . I n . > 1,(1  Til  PI 

I.  AI  L I'l  l.M  I I ¥l>l  I I , A I w , I I 
I c iiN  r I mir 

NV  A r ( V VT  - y MP  > / ( 1,0  . •(.!>  I r 1 ♦ I 
VM-  I Y MT  - VMII  I / ( NV  A- 1 I 
OH  I - I .NVA 

VIM' t = ( y Ml  1.'  ( V Mil  • y M*  I I _ I I I -V  PM  r I r 

VN'.ry  Mill  U -I  I •'.!  IO»I.O  . 

If|M()p(I.,'^|.li).0|  (,(.  TP  1^.. 

• ••.»ni-Aw  Hill  i/niiTAi  r.Pii)';  Aun  ANfuiTATr  III  mi.fifi;. 

(.All  III  T ( yl'  I l.ll  , VI'OT  . 1 I 
C A I I.  I'l  1 1 1 ( 'I . . y I 'I  I T , > J 
C A I 1.  n,  ,.11  I ' y I *1 1 T . Y>i  . , ( 
i.i' 

PH  CAll  I'll  (.(  I I ( y.-IlT  .YNS  , .I  1 
CAIl  i’li)T(n.,yopT,T| 
r A I L oi  U T ( yi'  I(,,i . ypiiT  1 
P2  C MM  I illif 

••*•*<’  (^  A'  o P.  ( II.  '■  - T.  » I c r I ’ foils  I. ONI. I Till  I ruPVFs  aim  oijawm. 

• THii.  »i-;r  (INI  y i.Fr.HiAii  mincatoi'  map  is  I'f'iivini  r). 

II  IIS1I  D.NI  ,ri  rAii  r ir  r ( I l_  AT  ,NS  .rtON.I  w ,■;<  tON.  SCI  AT  » 
yA«=-o,7‘, 

GO  TO  I A 

33  OI-INT  SS.yTilP 

6S  Igi-MAT(///<  STALi/.l/f  lf.|.  CI-.  PLOT  TOO  I AM;i  VImTICAHY.  VTnP=*. 

crr.y.*  ini  m s,  «a*  = to.o  iNCMrs.*» 

• ••••CAiL  INC'  I'l  oT  T(i  fi.osr  roi  puutinf 

7 CAl  L I’l  (Jf  ( r.  ,0.  ,QQ  >» 

STOP 
L NO 


OOOOG  A'fC 
OOOOSNOO 
corossi r 
ooonospo 
oooor.030 
or,ono<iAO 

onCflOQ*  0 

dOnorsAO 
OOOCOS  7C 
0 o 0 0 d s p n 

0 n 0 c d s <;  0 
00001 oco 

0000  I 0 I 0 

ooor I o?o 

0 O 0 '•  1 0 3 0 
0 0 0 p 1 0 A 0 

0 0 r 0 1 r ' 0 
00001 oco 

0000  10  70 

0000 1 OPO 

0 o 0 o 1 0 S 0 

0000 1 I 00 

oocoi  I I 0 

00001  I 20 

OdOOl I »0 

no 0" 1 I AO 
00001  1 so 

00  or,  1 I fO 

0000  I I 70 

00001  I PC 
OOOOl 1 oc 
OOOO I 200 
00001 210 
OOOOl 2?0 
OOOOl ? 30 
OOOOl 7A0 
OOOOl 2S0 

00001 ?eo 

OOOOl 270 


■.Hill-  OUT  I NA  r ir  T ( » L A T ,NS  .FtriM.I  W , SC  L ON  , !iC  L A T » 


this  sum  (M)T  I ni  is  r^  calculate  the  riCTiTinus  mepioians  and 

PAFALLfl  S SUPER IMPOSrO  ON  THE  MERCATOR, 


6C  001 2PP 
OOOO 1 200 
»*OO0Ol 300 
•OOCOI 310 
•OOCOI 320 
• 0000  1 3 30 
•OOOO I TAO 


CUM  MU  N XM.XMH.VMT,  yMH.YTOP.YROT.SSC.XRIOH.ISTITI 
OI  Ml  NS  I ON  NST  12  » ,Ey>r  ( A I , FPLUNC  IC.0  > 

INTI"  r, if  N/*Ni/.r/'r*/,NS,F*l,NST/*N*.»S*/.EMT/ir*,*l*»,»N», 
1 NTLGFP  SN.rLAr,  .Fwr  , I 0/0/,  KP/0/ 

Rf  AI,  RAO/.  01  7A5  12«)/.rCr./5  7 . 2 <J357  795/ , P I / 3 . 1 A 1 SRPftSA/ , MP  I / I 
Rsal  LrjN.LAT 

I F ( < XMT .CO. 1 OnOO.  I . ANP. ( XHP .FO.n,  11  KPsI 
I F < ( XMT  . r.f  .0000.1.  AND.  < XMR  ,LF  .-9900,  11  HP*? 

PRINT  H.FLAT .NS.FLUN.FY 

F ORMAT  ( |H*  ,r,X  . • F ICT  IT  I OUS  POLE  ■ • ,F  5.  2 , A 1 . 3X  ,F  ft  . 2 , A 1 // / 1 

IF  INS.  t 0.  Ml  C.n  TC  2 

FLAT=90. aflat 

GlJ  TlJ  3 

FL  A T = 00. -flat 

IF(n».Nf.F)  FL0N«360. -FLON 

I L ATP iPAO*  f LAT 

FLCNP  =PA0*FL0N 


•**».NO.  nr  POINTS  GENERATED  IN  OOTM  LONCITUOINAL  AND 

•••♦•latitudinal  directions 

NLON=( 360. /SCLON I 
NLAT  = ( 1 7S. /SCLAT I 
CH=SCLUNa60. *2. 

CPLONs-SCLON 
FLAG*0 

200  SCL ATaSCLAT/2, 

Lr*NLON 
LAxNLAT*  2»  I 

IFIFLAG.f O.I I CPLON»0. 

DO  10  1=1  ,LO 
CPLON=CPLON* SCLON 
FPLONI  I 1 sCPLON 

•••••F I C r I T I nUS  LATITUDE  CURVES  ARE  OPAWN  WHEN  FIAG»I 

IF (FLAG.EO.1 I GO  TO  203 
I F (FPLONI I I . GF , 1 HO . I GO  TO  38 
EMF  =F.I»T(  1 I 
GU  TO  39 


•••••••••OOCOI 3E0 

OOOOl =f 0 
OOCOI 370 
00001 3P0 
• S*/  00001  =90 

OOOOl 400 
, 5707Ae  3/00001  4 1 0 
OOCOI 420 
00001 430 
OOOOl 440 
OOCOI 450 
0000 1 4C0 
OOCOI 4 70 
00001 4P0 
OOCOI 490 
OOOOl EOO 
00001^10 
OOOOl T20 
00001 E30 
OOOOl SAO 
00001 S50 
OOOOl SCO 
00001 570 
0000! SPO 
00001 bSO 
OOOOl 6C0 
OOOOl 610 
00001620 
00001  ft =0 
00001 6A0 
00001650 
OOOOl fteo 
OOCOI 670 

00001 eeo 

OOOOl 690 
OOOOl 700 
00001710 
OOOOl 720 
OOOOl 730 
COOOI 740 
OOOOl 750 
00001760 


r>r>rw^nr»  r\r>n  nori  Ann  nnr. 


r ' r w M M 

1 PI  in  u » - w r . 

K,\)  T .1  ^ 

- I I'l  ONI  1 1 

^0^ 

If  (1  I’l  IINI  1 » . 

' .10.1  G(1  I n 205 

f »’i  UN  n ) - . • 1 

r - « w ’ ( 1 1 

; 1 ) T n * ■ > 

1 I'l  IIN  I I I 

f’O-i 

1 *1  M A 1 ( .4  ) 

r iM  '■!  r I ) 1 * t (' 

'It  1 1 - in  . 

' * 

i 1 . H *»•  »•  -M  • t»| 

c f'l  - 0, 

M 'N 

I M » t A(.  . f 0 . 1 1 

H R ► - 
K R 

R “ A 
Nw  - 0 

t)  1 .1  A 

KKR  --*i>».*  \ 

' I'l  A I - - SOL  A r 

I f If  1.  A(.  .f  IJ  . 1 ) 

f r = ► K 4 1 

C »M  A 7 -(  Pt  A 1 • 

i 1 A I 

C»M  All.  1 '•!  P| 

A = ( i '(  I 

4 T 

If  if  1 A'.,l  (J.  t ) 

M^C  *M  A T P 

4 -=CPl  ATP 

If  ( f 1.  A(,  • f (J  • 1 ) 

Hr  rPl  (If  l> 

•••••eOllAT  I('N‘.  11', M>  r ()*;  (.rNF'lAT  |N(.  nrTITIOUS  POINTS  IN  TMr  MrPCATOP 

Ai  =r<i'.(A)«ois(M  ATP)«ccr.  «rL  onp  > 

A,'  - S I fj(  A ) • >,  I N(  M.  (INP  I 

Alrf.lT  Afi((()*',  INIFLAIP  I*mS<  FLONP  ) 

AA-AI  -/.,  «AT 

HI  r I II'.  ( A I • CnSI  f 1 A TU  ( «S  INITLONP  I 
HJ  -M  N < A ) • CHS  ( n ONP  » 

rll  - ( I II  Arj  ( fi  1 • S I HI  rt  A I 0 » «S  IN  irtONP  I 
n,i  = (i  1 4 II  4ii  1 
c = ‘.HI-  I < A A*  •?  4HM«  »2  I 

0 = I - rO‘.  ( A I *S  I N(  f I,  ATP  I ♦ COT  AMR  I •rosin  ATP  M 
CAT  ArAN?(r,o» 

LON- A I ANP ( CH , A A I 

••• **0 ALrui  A Tf  THf  latituof  nr  point  o in  MFRCATOP  systfm 

1 I (I  A T ,(.t  , 4-0  II  or  TO  TS 
l,A  T iHf  I -la  t 

AL  A T tL  A I • AO.  ♦I'OC. 

SNiNST  I I t 
C4O  T . I I A 

35  LATilat-HPI 

AL  AT  -r-L  A T«r  0 , *l>l.  G 
SN^NST I 2 I 

36  OLA  T =LA  T «(if  (. 

• • • . • CAL  I III  A T r TMF  LONOITUOF  OF  POINT  O IN  MFPCATOP  SYSTFM 

IF  I I (IN. GI  .0.  I GO  TO  17 
Cfc  =f  wT  I 2 I 

HL(  N- Acsil  oimoro 

ALON  = I OHOr  . -HI.ONA60  . 

GO  TJ  I .• 

3 7 CW=(  WT(  I I 

HLCHU  (IN»f'FO 
ALON^Ml  nN»A.n.-I0«0O, 

1 I I I J.  I.I  . 1 I on  II)  IT 
ALf  M rALf  N 
ALA',I  =ALAT 
13  ALijsr -ALUS I 
AL  AS  2 dalasi 
ALPS! lALCN 
ALAS  1 ^ALAT 
FPT  =XMT-ALrS2 
F P 1 - * MM- AL  ns  2 
IF  I K.I.r  . I I r.I)  TO  16 


0 0 00 1 A 70 

0000  I 7F0 
00001 700 

oconi  1-00 

00001  PI  0 
00001 0?C 
00001  P TO 
r 0 0 0 1 p 4 0 
0 0 0 0 1 ' 0 
0000  1 417  0 
0 0 0 01  •’  7 r 
0 0 001  P PQ 
OC'OOii  p<;o 
00  001  '!00 

00  00 1 010 
00001 9?0 
OOOOl S 30 
00001 940 

00001  9*-0 

0000 1 OF 0 
OOOOl 970 
0 0 0 0 1 o n 0 
ooor 1 990 
00on;o00 
00002010 
OOOOJOJO 
000020 10 
0 0 0 0 ? C 4 0 
00c 02050 

00  002  OF  0 

00 002070 
0 0 0 0 2 0 A 0 
0 0 0 0 ; c-  9 0 
OOOf  2 I 00 
0 0 00 2 I 10 
000021 20 
cor  021  30 
000C2I 40 

00002  I '0 

00002  1 eo 

COOOi 1 70 
00002  1 AO 
00Ci'21  90 
000022  00 
00002210 
0O0022P0 
00002230 
OOr>0224C 
00002250 
00002270 
00002270 
onoorpeo 
00002290 
''’''O23C0 

0000;  . i 0 

00C02320 
00002  3 30 
C0002340 
00002350 
00002370 
00002  3 7u 
000023PO 
00002 1 90 
on  CO 14  00 
0000241  0 
00002*20 
00002*10 
00002**0 
00002*50 
00002*60 
00002*70 
00002*80 
00002*90 
00002*.00 


• •••*CMfCF  IF  POINT  IS  OUTSIUE  TMK  MPPCATflR  MAP 

iri  Al  AT. (,T. YMT.OF.ALAT.lt. YMIll  GO  TO  21 
IF  I A 1 ON.r.T  .XMT  .OB.  ALON.LT  . KMOl  GO  TO  21 
1 F I J.  r O.  I ) f.n  TO  5 
l»  I K .1.1  .1  1 GO  TO  16 

• ••••r  XTPAI’I  I ATI  Till  COI'Vr  TO  T»4K  FOGF  if  data  point  is  roUNO 

• •*  ••OUT  s I nr  Till  MAI’ 

• ••••SPICIAL  CAS'-,  will  N PIOTTING  IS  HANOLFII  FROM  OW  TO  OF  OR  OR  TO 

IF  I KP  .1  0.0  1 I4I'  TO  71 

|l  1X1*.  f,(  .1  .AIIIi.NW.f  I). I I FtO  TO  7 2 


00002510 

00002520 

00002510 

000025*0 
00002550 
0000257 0 
00002570 
00002  5 80 
00002590 
00002700 

mow  00002610 

000027  20 
00002730 
000027  *0 


i 


0. 


I • ( At  ll'..'*  .<•(  . 0 . , * Kf  . M I !«,?  .1  r . KMT  .AND.  At  I tN  . I » . )l  MT  . A NO  . Al  riM.  r,|-  , Q.  » GO 
(.III  II 

If  < A I I r,.'  , I { . 0 . . A Ml.  Al  1 I*-..-’  ,(,r  ,KMM.  AMO.  Al  ON  .1  ( ,0  . , AND  . A|  0(1  , «M0  »Gn 

T I ) T } 

M AIIM  A|  ( .1  ( .cm.aNO.  AOSI  Al  ONI  .1  r .0(1  GU  To  71 

If  I Al  n<- ,>.  (,r  . y MM  . a no  . A|  ON  . l.  1 . 0 . . ANO  . A|  nr.  ? .1  I .0  . I A|  rvir’-'KiAT 

If  ( ► o . f o . . AMU.  Al  t.'w’  .1  r . X MT  . Ann.  Ai.rM  .1  t . o . I ai  os?  -*mo 


(.1 ' 

T • 1 

n 

7? 

If 

( Ai 

u s . r 0 . 

XMT  . 

IIP  . Al  (1 

S?  ,7  C . X MM  1 NW=0 

|7 

( Ai 

r'*./.  f ‘1. 

XM  T > 

(.11  Tn 

77 

I f 

( Ai 

H M M 1 

Al  ns? 

- X M T 

G'  1 

T • 1 

7| 

7 7 

Al 

1 

n M 

7 1 

If 

( Ai 

A ■ .<.f  • 

YMT  1 

C At  L 

7tiGFPI  1 Al 

AS?.Al  OS2.AI  AT  ,A|  fIN.PI  A S , PI  f)S  . 1 1 

II 

i Ai 

.1  r . 

X Mil  ) 

f Al  1 

1 fK.rpi  I Al 

AS2  , Al  OS?  , Al  AT  , Al  nN.PLAS  .IN  ns.  2 1 

If 

( Ai 

A'-..' . 1 r . 

Y M ) 

CAI  t. 

1 OGI  PI  < Al 

AS2.AI  nS2.  Al  AT,  ALPN, PI  AS. PI  iIS,  1 1 

17 

( Ai 

1 's.'  • (.1  • 

X M 1 ) 

CALI 

rDG7P1  (Al 

AS?.Al  OS2.AI  AT  .ALON.I’LAS.PI  PS.  21 

X T 

► ♦ 1 

Al 

« fJ| 

- ( fi 

X M i U 

•SSC 

AL  A I 1 - ( YM(  I.  ( I ( ,\s  l-yncT  I • ssc 

C Al  L I’l  01  ( At  IINl  . Al  ATI  , 1 I 

r,n  TO 


ooor  ? ( 'c 
Of’OO?!  to 
Ooon?t  70 
0000?  t AO 
O O O n ? ( (.0 
0 0 00?  700 

oooo;>  71  0 

0000?7?0 
0000  ? 7 '0 

ooour  7AC 
onoo  ? 7<-  0 
0000?  77  0 
0000?  7 70 
0000?7('0 
0000?700 
ooof  ? f'CO 
0000? CIO 
0000?c?0 
OOOOPH  10 
oooo?oAo 

0000?PAO 

0 0 00?  0 7,0 


c 

C*«**«ri<FrK  |l  fflMT  IS  rniTSIOF  1(17  MFPOATnP  MAI> 

C 

|(S  If  I Al  At  ,(,T  .Vdl  .llG.Al  AT  .1  T . vxni  GO  TO  ?l 
I r < Al  ON.  gT  . « MI  . OP.  Al  TN.l  T.KMHI  GO  TO  ?l 
C 

c • Sf'f  C I Al  (ASf  MfOCATHC  plot  rP(1M  O W TO  07 

C 

17  ( Kt>.  f O.  ?.  AMf'.r  r«T  ,LF  .cm  . ami  .Al.ON.r.7  .XMH.  ANO.  ALON.t  T . 0 , | Al  ON=  *MT 
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